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LANDING SITE DISPERSIONS OE AN UNCONTROLLED 


LANDER ON MARS 
By John S. White 
Ames Research Center 


SUMMARY 


A study was made to determine the accuracy with which a lander, descend- 
ing from a bus orbit, can reach a specified landing site on Mars without using 
a range control system. This report considers both a lifting (l/d = 0.5) and 
a nonlifting lander - each with two different values of mass. Several differ- 
ent entry angles were considered, along with two different procedures for the 
deorbit maneuver; the use of the minimum deorbit velocity increment and the 
use of a constant deorbit velocity increment. 

The error sources considered are variations in the Martian atmosphere, 
variations in the actual lift of the vehicle, and errors in performing the 
deorbit maneuver. It is concluded that if minimum deorbit velocity increment 
is used with entry angles of -20° to -24°, the landing site accuracy will be 
on the order of 1°, which is of the same order as the present knowledge of the 
location of the landing site. Thus a range control system need not be 
included in the lander. 


INTRODUCTION 


One plan for the unmanned exploration of Mars is to use an orbiter bus to 
carry a lander into an elliptical orbit. After several orbits, during which 
the knowledge of the orbit is improved, the lander is detached and given a 
small AV to put it in a descent trajectory. The vehicle then enters the 
atmosphere and experiences some atmospheric braking. Just before touchdown a 
soft landing device, either a parachute or a retrorocket, or both, is used. 

The landing trajectory will be chosen so that the vehicle will touch down at 
some specific location on the planet. 

Comparisons of landing vehicles with and without lift have shown (refs. 1 
and 2) that if a lifting vehicle is used, the drag-to -weight ratio can be 
reduced, or equivalently, the landed payload can be increased. (Some consid- 
erations of entry requirements for lifting and nonlifting vehicles are given 
in refs. 2-4.) 

In general, we can expect that the landing vehicle will not touch down at 
exactly the desired point, but instead the touchdown point will have some dis- 
persion. If this dispersion is sufficiently small, then the landing will be 
satisfactory. If, however, the dispersion is too large, then some form of 
control must be used during the descent trajectory. This control might be 


applied either immediately after the application of the deorbit AV or else 
during atmospheric flight. The purpose of this report is to determine the 
size of the dispersion when no control is used, and from this data, to deter- 
mine whether a control system is required. It is desirable, of course, to 
avoid the use of a control system, since without one the vehicle would be 
simpler and, therefore, more reliable. 

The causes of landing errors considered are: variations in the Martian 

atmosphere, errors in making the velocity maneuver to enter the descent tra- 
jectory, and uncertainties in the lifting characteristics of the vehicle. 
Pragluski (ref. 3) did work similar to that reported here, but he did not 
carry the trajectory through the atmosphere, and used only a single entry 
velocity. Avco (ref. 4) has also done similar work, but restricted it to a 
nonlifting vehicle, a single entry angle, and a restricted set of entry ranges. 
This report considers the more general case of a lifting or nonlifting vehicle, 
with two different values of mass, entering the planetary atmosphere anywhere 
in the plane of the initial orbit, at several different entry angles. The ini- 
tial orbit is assumed to be equatorial and posigrade. The effects of other 
orbit inclinations are considered in appendix A. The deorbit AV used is 
either the minimum value for entry at a particular range, or a fixed value of 
sufficient size to allow entry at any point. For both cases, the expected 
downrange error at impact with the surface is determined. The report presents 
results for the planar problem, and some theoretical considerations of out -of - 
plane motion are included in appendix B. The effects of Martian winds and of 
motion while the soft landing system is in operation are not considered. 

In order to determine whether a particular landing site error is accept- 
able, the error must be compared with the available knowledge of the Martian 
surface. The knowledge of the location of desirable landing sites will be 
based on astronomical photographs, and on photographs taken from spacecraft 
such as Mariner IV or the orbiter bus itself, prior to detaching the lander. 
Slipher (ref. 5> P* 59) > in discussing the quality of earth -based photographs 
of Mars, says "The finest linear markings found on our best images measure 
around O'.l to 0l ! 5 of arc, but their true width must be greater because irradi- 
ation will tend to narrow them down considerably. The smallest circular 
objects are about three times this diameter for dark spots, and even greater 
in case of white spots.” These remarks refer to resolution, that is, the abil- 
ity to distinguish between neighboring objects. The ability to specify the 
location of a particular feature is not that good. Thus, the location of 
desirable landing sites is known no better than about 0V2 of arc on Mars as 
seen from the Earth during close oppositions. For comparison this is equiva- 
lent to about 50 km over the surface of Mars, or about 1° of central angle. 

The Mariner IV photographs have much greater resolution, that of the best of 
them being about 1 km (ref. 6). However, the attitude control system of 
Mariner IV had an accuracy of about 0.5°* Since the closest approach was 
about 6000 km, the uncertainty in the location of the picture is again about 
50 km. If the picture is taken from the orbiter itself, then some improvement 
in the knowledge of the location of the picture .can be achieved. If we assume 
the same attitude control accuracy, and a 1000 X 5000 km altitude elliptical 
orbit, then at perigee the uncertainty of the picture location will be about 
10 km, whereas at apogee it will still be about 50 km. Thus, it seems reason- 
able that for many situations a landing error of 50 km, or 1°, would not be 
considered large, and even 2° error may be acceptable. 
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NOTATION 


f vehicle acceleration per unit mass , km/sec 2 

fy lateral component of aerodynamic force per unit mass, km/sec 2 

i inclination of orbit 

— lift -to -drag ratio of lander 

ballistic coefficient, kg/m 2 

CpA 

r radial distance from Mars, km 

velocity on descent trajectory at intersection with bus orbit, km/sec 
V e entry velocity with respect to inertial space, km/sec 

Y-^ out -of -plane component of landing site, km 

7 e flight-path angle at entry with respect to inertial space, positive 

upwards , deg 

&( ) variation in ( ) 

AV magnitude of velocity vector required to transfer from bus orbit to 

descent trajectory, km/sec 

AVy out-of-plane component of transfer velocity vector, km/sec 
r) angle between AV and local horizontal, deg 

0 range angle, deg 

0 e range angle of lander at entry, measured in direction of motion from 

periapse of bus orbit, deg 

0^ range angle of landing site, measured in direction of motion from 

periapse of bus orbit, deg 

£ out-of -plane component of angle between AV and local vertical, deg 

cp vehicle roll angle, deg 

out -of -plane component of landing site, deg 
( ) sa three sigma dispersion of a statistical quantity 
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ANALYSIS OF RANGE ANGLE ERRORS 


The problem is to determine the expected range angle error on landing, as 
a function of the three principal error sources: errors in making the veloc- 

ity maneuver, variations in the atmosphere, and variations in the vehicle f s 
lift-to-drag ratio. Here, range angle is defined as the angular distance 
traveled over the surface of the planet in the plane of motion. 

Navigation errors (i.e., initial condition uncertainties) at the start of 
the entry are not considered. After several orbits of Mars, the expected 
value of these errors will be on the order of kilometers in position and 
meters per second in velocity. The position error propagates to the surface 
nearly 1:1, and is negligible compared to the other errors. The velocity 
error is smaller than the error in making the velocity correction, and 
therefore can either be ignored or combined with the velocity correction. 


Let us arbitrarily break the trajectory into two parts, the descent tra- 
jectory, immediately after the velocity maneuver and continuing down to the 
entry altitude, and the atmospheric entry, from entry to touchdown. This divi- 
sion of the trajectory lessens the 
number of parameters to be con- 
sidered in each part, since the 
first part is independent of the 
vehicle and atmosphere, and the 
second part is independent of the 
bus orbit and the range angle at 
entry. The sketch shows the geome- 
try involved and defines a number of 
parameters used in this analysis. 

From the first part of the tra- 
jectory, one obtains entry condi- 
tions of V e , 7 e , and 0 e as a func- 
tion of the velocity maneuver. To 
a first-order approximation, devia- 
tions in these entry conditions are 
given by the expressions: 
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(i) 


Geometry of entry maneuver. 


where 5AV and brj are, respectively, 
the magnitude and direction errors 
in making the velocity maneuver. 
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From the second part of the trajectory, for a given atmosphere and 
vehicle, one can obtain range angle as a function of entry conditions. Using 
the first-order approximation as before, we have 




+ SV a + 6 ?e 


av. 




( 2 ) 


Since the altitude of entry is also an entry condition, we might expect 
this to appear in equation (2); and if the initial point of the trajectory is 
specified as a particular time, this would be true. However, we have speci- 
fied the end point of the descent trajectory as that point at which the tra- 
jectory passes through the entry altitude, regardless of whether the vehicle 
is in a perturbed trajectory or not. As a result, of course, the partials of 
equation (2) are not dependent on altitude, but must be calculated with this 
particular starting condition. 

Equations (l) and (2) combined give the range angle error due to errors 
in making the velocity corrections. 

The errors in range angle due to atmospheric variations were obtained by 
numerical integration of the equations of motions, starting with a fixed set 
of entry conditions and using different atmospheres. A comparison of Qj^ for 
the different atmospheres then determines the deviation due to atmosphere, 

^^atmos • 

Finally, the first-order error due to lift -to -drag variations was 
determined as 


* _ 



( 3 ) 


The error sources, 5AV, 6r], 5 (l/d), and the atmospheric variations are 
independent. The overall 3cr expected landing site error can then be 
obtained by adding the individual 3<? errors on a root-sum-square basis, 
giving 


Y^L ^e) &AV I 2 

_W e bAV by e bAV be e bAVj 30 _ 


a / 2 


(&e L ) 


total, 


30 


, dv e $>y e ae L be e \ 
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The quantity 50^ was obtained from the expected extremes of the possible 

a ci Dios 

atmospheres, and it was felt that it could be considered as a 3c number; 
thus it is of the appropriate size to combine with the other error sources. 


COMPUTER SIMULATION 


In order to determine the numbers to be used in equation (4), two simu- 
lation programs were written for use on an IBM 7094 computer, one for each 
portion of the trajectory. It was assumed that the descent trajectory and 
entry were posigrade equatorial orbits about Mars, so that only two-dimensional 
programs need be generated. A brief study of out-of -plane errors and of the 
effects of orbit inclination are included in the appendices. 


Descent Trajectory 


The first of these programs considered the descent trajectory, between 
the velocity maneuver and entry into the Martian atmosphere. This program was 
used to obtain the descent trajectories, and the partial derivatives of condi- 
tions at entry, V e , / e , 0 e , with respect to the magnitude and direction of the 
corrective maneuvers, AV and 3 (eq. (l)). The procedure used in the program 
was as follows: First a bus orbit was selected. The orbit altitude was 

5000 km at apoapsis and 1000 km periapsis. With the bus orbit selected, a 
descent trajectory was determined which used the minimum AV, 1 and which 
arrived at an altitude 2 of 244 km with specified values of entry angle, y e , 
and entry location, 0 e . This calculation was repeated for -l80° < 0 e < l80°, 
to allow considerations of landing at any point in the plane of the orbit. 

Zero range angle was chosen under the periapsis of the bus orbit. The entry 
angle was also varied from -24° to -l6°. The entry velocity, V e , was deter- 
mined for each of these trajectories, and plots of A Vm-i n and V e vs. 0 e 
are presented in figures 1 and 2. 


For operational reasons, it may be desirable to use a fixed value of AV, 
which will still allow landing at any location. The value chosen, AV^^^, 
would be the maximum of the AV m j_ n shown on figure 1. The particular values 
used are given in the following table. 


7e 

^min max 


-16 

0.94 


-20 

1.07 


-24° 

1.24 km/sec 


Thus, a trajectory was determined for each 0 e and y e which used the speci- 
fied AV m in max * The corresponding V e was also determined, and is plotted in 
figure 3* 

1 It is interesting to note that AVm-f n corresponds to the minimum V e . 

2 The entry altitude was chosen as 244 km since this is just beyond the 
outer reaches of the highest atmosphere considered. The descent trajectory 
was stopped and the atmosphere entry was started as the vehicle passed through 
this altitude. 
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Having determined the trajectories to be considered, we next compute the 
partials of 0 e , / e , and V e with respect to AV and t) . This is done by using 
the two sets of trajectories just determined and perturbing AV and t] slightly. 
The partials are then computed from the perturbed trajectories. Standard 
elliptical equations were used throughout this program. 


Atmospheric Entry 

The second program was used to determine the effects subsequent to entry. 
In this case, a single trajectory was integrated from the desired entry condi- 
tions (v e , 7 e ) to touchdown, where V e and / e were obtained from the trajec- 
tories determined in the previous section. The vehicle is assumed to have a 
constant ballistic coefficient and a constant lift-to-drag ratio throughout 
the entry, with the lift applied in the upward direction. Entries were made 
using two different values of each of these parameters. The ballistic coeffi- 
cients used were m/CpA = 62.9 and 23.6 kg/m 2 (0. 35 and 0.15 slug/sq ft) and 
these cases are referred to as the high and low mass cases, respectively. 

The values of l/d used were 0 and 0.5. The equations used in determining 
the entry trajectory are those of reference 7 * 

Thus a family of trajectories was obtained, from which range angle was 
plotted versus entry flight angle, initial velocity, and l/d. The slopes of 
these curves were then used to determine the partials of range angle with 
respect to velocity, flight -path angle, and l/d for a particular atmosphere. 
In addition, entries into different atmospheres, with all other parameters the 
same, were compared to obtain the dispersion due to uncertainty about the 
atmosphere. The atmospheres used in this program are assumed to be exponential 
above some transition altitude, and convective below. They were based on the 
JPL atmospheres VM3, 4, J , and 8 (see ref. 2). It was also assumed that the 
atmosphere was rotating with the planet and that there was no wind relative to 
the surface. 


Landing Conditions 


It is of interest to consider the landing conditions that will exist for 
the variety of vehicle parameters and entries under consideration. In partic- 
ular, we would like to be sure that the vehicle becomes subsonic at a suffi- 
ciently high altitude, so that the soft-landing system will operate properly. 

Accordingly, the altitude at which the Mach number was 0.9 was determined 
for all entry trajectories and vehicles, and for the two extreme atmospheres. 

It was found that this altitude did not vary significantly with entry angle or 
velocity, but did vary with mass, l/d, and atmosphere as shown in the following 
table: 
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Mass 

Atmosphere 

Altitude, km 

l/d = 0 

l/d = 0.5 

Low 

VM8 

5 

11 

High 

VM8 

-4 

3 

Low 

VM3 

22 

27 

High 

VM3 

7 

Ik 


From this it can he seen that the high mass vehicle with 0 lift is not 
practical, since on entering the VM8 atmosphere it will impact before slowing 
down sufficiently. However, in order to present a complete picture, the 
results of this case will be presented along with the others. 


RESULTS AND DISCUSSION 


The results are given in two parts. First the effects of the individual 
error sources on landing site dispersion will be presented. Secondly, these 
errors will be combined to give a picture of the overall dispersions. 

The four error sources under consideration are, as mentioned earlier, the 
error in the magnitude of the deorbit impulse, 6AV, the error in the direction 
of the deorbit impulse 5r], the error in the lift characteristic of the 
vehicle S(l/d), and the error in the knowledge of the atmosphere at Mars. 

The particular values used here for these error sources are 


5AV scr = 0.75 percent 

5 ^sa = 1 ‘5° 

5(L/D) sa = 0.05 


while the atmosphere error represented the extreme error caused by the atmos- 
phere uncertainty. These numbers are 3^ values of the expected error. The 
first two were obtained from reference 5 (p* 91) j and the third number is 
10 percent of the maximum L/d of 0«5» The first set of results, showing the 
effects on the landing site error, SDp, caused by each of the error sources 
individually, are plotted in figures 4(a) -(h). The partials plotted in these 
figures are for entries into a VM8 atmosphere. Partials were also obtained 
for entries into a VM3 atmosphere, but the curves are quite similar, and are 
not presented. The dispersion due to atmosphere, however, is the extreme over 
all four atmospheres. 

Several trends are apparent from the curves of figure 4. When A Vmi n 
is used (figs. 4(a) -(d)) with l/d = 0, the miss caused by 5 (l/d) is negli- 
gible. At L/d = 0.5> however, S(l/d) can cause significant errors. This 
increase in miss due to S(l/d) with L/d is independent of the vehicle mass. 
Also, at L/D = 0.5, the effect of atmosphere is negligible for 7 = -24° and 
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low mass. It increases somewhat with mass, but much more so with decreasing 
entry angle. This increase with decreasing entry angle is general for all of 
the partials, and at the shallower angles any one of the four errors may be 
the larges t, depending on the value of 0 e . 

For AV-m-f n m32 r (figs. 4(e) -(h)) the same general remarks hold, with the 
exception of the errors due to Arj. These errors are much larger than any 
other in the region from 0 e = to = “30 j which is the region where 

AVmin is considerably greater than A V Trn < n . 

The overall results are shown in figures 5(a) -(h) where the 3a value of 
the total landing site error, 50L^. o ^ a j_^ is plotted against entry range# The 
curves for atmosphere VMS are the root sum square of the curves in figure 4* 
The curves for atmosphere VM3 are obtained from data similar to that of fig- 
ure 4, in which the partials were computed for entry into the VM3 atmosphere, 
whereas the dispersion due to atmosphere variations was taken exactly from 
figure 4. Thus, the curves show that there is little difference in the dis- 
persion when partials are calculated for the VM8 or the VM3 atmospheres. 

From these curves one can also conclude that: the error is considerably 

larger for the shallower entry angles, it increases slightly with L/d at the 
steeper entry angles, and more so at the shallower entry angles, and is nearly 
independent of mass. 

If A V-mi n and the steep entry angle j - -24° are used, the 3a entry 
error will be less than 1.5° for any entry range, and less than 1° for some 
entry ranges. Using AV-mi however, may lead to a restriction in the 
choice of entry ranges, since the error is less than 1.5° for only these entry 
ranges between 0 e = -3O 0 and 0 e = -120°. If it is desired to remove this 
restriction, some control system will be needed. Since the only large error 
component in this region is that due to br\ (figs. 4(e) -(h)), it is sufficient 
to control only this error source. Thus, a system which measured rj, and 
subsequently used a vernier thrust to correct for errors in rj would be 
sufficient. 


AV 


It is interesting that in the region of 0 e = 0, the error associated with 
is less than the error associated with AVn,^. The reason for this 


-q _l_l^ UiiUllJ. -L V-/J. UiU UUiW V1J. UJ.J, LJ V n * 

max 

is that | ^ > 0j J /^ >T l | is smaller for term is the dominant error 

source in that region. 


Similar data obtained for other orbits of different eccentricity indicate 
that the same general trends occur, and that similar conclusions can be drawn. 


CONCLUSION 


It has been shown that a landing vehicle entering Mars from a bus orbit 
with entry angles of -20° to -24° will land within about 1.5° of the desired 
landing site without the use of a range control system after deorbit. It has 
also been shown that at present our knowledge of the location of the desired 
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landing site is approximately 1° to 2°, so that a landing site accuracy of 
that size appears quite reasonable. 

If the minimum AV is used for the deorbit maneuver, then any location 
in the orbital plane can be reached with that accuracy. If, however, a fixed 
AV, equal to the maximum of the possible minimum AV is used, then the land- 
ing sites which can be reached with an accuracy of 1.5° are limited. The land- 
ing sites available, for this condition, are those between 30 ° before perigee 
of the bus orbit and 120° after perigee. If it is desired to use the 
AV mi n max , and still require that any landing site be achievable, some correc- 
tive scheme must be used. An adequate scheme would be one that would reduce 
the error in the direction of the velocity correction by a factor of about 6, 
from 1.5° to 0.25°, 3a* The other error sources are sufficiently small that 
they need not be corrected. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035 * March 4, 1968 
125-17-05-13-00-21 


APPENDIX A 


EFFECT OF ORBIT INCLINATION 


The main study has assumed that the vehicle will be entering the Martian 
atmosphere from an equatorial posigrade orbit. The effect of different incli- 
nation orbits will be discussed here. There are two different effects to 
consider. 

First, as the inclination varies, the entry velocity with respect to the 
atmosphere will vary, causing variations in the flight time, and a resultant 
in-plane error of the landing site. The two extremes are posigrade and retro- 
grade orbits, and the dispersion due to different atmospheres was studied for 
these two cases. For entry angles of -20°, this dispersion was greater for 
the retrograde orbits by about 0.2°. The change in dispersion was less for 
steeper entry angles. This is only a modest change so that the general results 
shown in the text apply for any inclination angle. 

Second, any variation in time of entry will cause a landing site error 
due to the rotation of the planet. For an equatorial orbit, either posigrade 
or retrograde, this effect has been considered, since the landing site has 
been computed in Martian-fixed coordinates. For inclined orbits, however, 
there will be an out-of -plane component, with a maximum for i = 90°. The 
maximum deviation in time of entry due to atmospheric variations is about 
150 sec, with a corresponding maximum lateral miss of about 0.6°. 

Thus, the effect of orbit inclination on landing site error is relatively 
small. 
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APPENDIX B 


OUT -OF -PLANE ERRORS 


The main study considers only the affects of in-plane errors. However, 
a brief study of out-of -plane errors was also made. Two possible out-of -plane 
sources were considered. The first was an error in the direction of the veloc- 
ity maneuver, and the second was an error in the roll attitude of the vehicle 
during the entry. 


Velocity Manuever Error 

If the velocity maneuver has an out-of-plane error, the plane of the 
actual descent orbit will have an angle of inclination i = sin-i(AVy/v^) , 
where AVy is the out-of-plane component of the velocity maneuver, and 
is the velocity of the descent orbit just subsequent to the maneuver. We then 
further assume that the vehicle remains in the new plane. The out-of-plane 
distance at the landing site due to this inclination is 


Yp = r sin i sin(0p - 0 e ) 


where ($p - G e ) is the angular distance traveled along the descent trajectory 
to touchdown, and r is the radius of the descent trajectory at the point 
where it passes over the landing site. We are interested in obtaining an 
approximation of the maximum value of Yp. In so doing, we should use a large 
value for r, such as the entry radius, r = 3^72 km. We can write sin i as 
AV/V d sin 5|, where S£ is the out-of -plane component of the directional 
error br\; and we can use the maximum for sin(0p - 0 e ), namely, 
sin(0 L - 0 e ) = 1.0. 

With these assumptions and for = 6 t] = 1.5° 


= 3472 vj (0 - 026) 

Converting this to an angular displacement over the Martian surface, in 
degrees, we have 


t 


^max 


3472 (AV/V d )( 0.026) 
3332 


57-3 « 1.5 


A V 

V., 


For trajectories with A V ^-; n , the maximum value of AV/v d , considering 
all entry ranges and angles, was AV/V d = 0.5. For trajectories with 
A Vmi n max this maximum value was AV/V d = 0.66. The resulting maximum 
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out -of -plane error due to the velocity is thus less than 0.75° with A V^j n , 
and less than 1° with A 

ULLL 1 msLX. 

Roll Attitude Error 

If , however, the out-of -plane error is caused by a constant vehicle roll 
error cp, there will be a lateral acceleration, f y’ proportional to the 
vehicle acceleration f; that is 


f = f sin cp 

y Y 


Here it is assumed that f is the total vehicle acceleration, and that 
the vehicle roll axis is perpendicular to the total acceleration vector. The 
lateral displacement at landing will then be 


, landing 


y l=, 


f sin cp dt 2 = sin cp 


, landing 


f at* 


entry 


entry 


This integration was performed for two different trajectories. Both had 
a shallow entry angle, 7 = -12°, resulting in very long flight time, and dif- 
fered only in the atmosphere used. Thus the results give the maximum expected 
lateral range. The results were that 


i|r T = 0.6° for VM3 atmosphere 
mnax 

\|r T = 0.42° for VM8 atmosphere 
J-'max 


If the results for the two error sources are combined on an rss basis. 


the maximum lateral range will be less than 0.9 with AV, 


1.2 with AVi 


mm; 


and less than 


^“max* 
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Figure 1. - Minimum AV as a function of entry- 
range and flight -path angle. 
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Figure 2.- Entry velocity associated with AV^n 
as a function of entry range and flight-path 
angle. 
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Figure 3* “ Entry velocity associated with the maximum as a function of entry range 

and flight -path angle. 
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Figure 5-- Concluded. 
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